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ABSTRACT 

Galactic transport models for cosmic rays involve the diffusive motion of these particles in the interstellar medium. Due to the large- 
scale structured galactic magnetic field this diffusion is anisotropic with respect to the local field direction. We included this transport 
effect along with continuous loss processes in a quantitative model of galactic propagation for cosmic ray protons which is based on 
stochastic differential equations. We calculated energy spectra at different positions along the Sun's galactic orbit and compared them 
to the isotropic diffusion case. The results show that a larger amplitude of variation as well as different spectral shapes are obtained in 
the introduced anisotropic diffusion scenario and emphasize the need for accurate galactic magnetic field models. 

Key words, cosmic rays, astroparticle physics, diffusion, ISM: magnetic fields, Methods: Numerical 
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1. Introduction 

Modelling the cosmic ray (CR) transport in the galaxy is a fun- 
damental topic in high energy astrophysics. Such studies are of 
great importance in the analysis of CR origin like supernovae 
or pulsar winds. The characteristics of the models depend on 
the properties of the interstellar medium (ISM) through which 
the particles travel and can thus give some insight into its fun- 
damental constitution. For example, the galactic magnetic field 
and its turbulent component are connected to the transport pa- 
rameters in such models. A reference on general CR propertie s 
can be found, e .g., in the tex t book s by Berezins kii et alj (|l990), 
iGaisserl (fl99(il) . ISchlickeiserl d2002l) or lStanevI (l2004l) ~ 

The basic propagation process of CR in the ISM is the 
diffusive motion of the particles due to scattering at mag- 
netic field fluctuations. From numerous studies in heliospheric 
physics, it is well known that the diffusive transport of ener- 
getic particles cannot be described by a scalar diffusion coeffi- 
cient but requires a diffusion tensor which takes into accou nt that 
parallel and perpend icular diffusion are different (e.g. Jokipii 
ll966r.lPotgieterll2011h . In galactic propagation studies, however, 
anisotropic diffusion of CRs has been investigated only for ba- 
sic magnetic field configurations o f partly l ocalized applicabil- 
ity, se e, e.g.. IChuvilgin & Ptuskinl d 1993b . iBreitschw erdt et al. 
(2002), Snodi n et alj d2006l) andreferences therein. While the 
latter authors were interested in the consequences of anisot ropic 
diffusion for ene rgy equipartition, lHanasz & Leschl d2003l) and 
iRvu et all d2003l) analyzed its importance for the Parker instabil- 
ity- 
More recently, Han asz et al.l d2009l) found that anisotropic 
diffusion is an essential requirement for the CR driven mag- 
netic dynamo action in galaxies. Their assessment of the dif- 
fu sion anisotropy is ba s ed on the full-orbit analysis performed 
bv lGiacalone & Jokipiil d 19991) who found the perpendicular dif- 
fusion to be significantly lower than the parallel one in a broad 
energy range and for both, isotropic and composite turbulence. 
Moreover, a recent derivation of the perpendicular diffusion co- 
efficient for galactic propagation, using the enhanced nonlin- 



ear guiding center theory and a Goldr e ich-Sr idhar turbulence 
model was performed by Shalch fet al.l d2010l) and resulted in 
ratios between the parallel and perpendicular diffusion coef- 
ficient which where much lower than unity as well, namely 
k±/k\\ ~ 10~ 4 - 10 _1 , depending on particle rigidity. An analysis 
for different turbulence spectra i n the context of super nova rem- 
nant shock acceleration of CRs (Marcowith et al.ll2006|) or their 
transport in chaotic magnetic fields dCasse et al.ll2002l) yielded 
similar values for this ratio. 

Many popular models for galactic CR transport, however, 
inclu de only a single dif fusion coeffic ient, like t he GALPROP 
code dStronget al.ll2010t) . Although Istrong et all d2007l) princi- 
pally acknowledge that anisotropic diffusion is of importance, 
they argue that due to large-scale fluctuations in the magnetic 
field on scales of the order of lOOpc, the global diffusion will be 
spatially isotropic. Observations of the galactic magnetic field 
indicate though that the field has a large-scale ordering with an 
regular field strength with ab out the same magnitude as the tur- 
bulent component (see, e.g.. lFerrierel200ll) . A similar indication 
is given by observa tions of external spiral galaxies (Beck l20TTt 
iFletcher et al1l201 ll) . which show a global magnetic field struc- 
ture aligned to the spiral arm pattern. Therefore, it must be con- 
cluded that anisotropic diffusion can have an important effect in 
galactic CR propagation. 

Besides from its fundamental astrophysical relevance, the 
spatial distribution of CR flux in the g alaxy is also of interest 
in the context of long-term climatologv. lShavivl d2002) proposed 
a CR-climate connection on the timescale of 10 8 years due to the 
transit of the solar system through the galactic spiral arms during 
its orbit around the galactic center. The argument assumes that 
the low-altitude cloud coverage increases due to an increased 
formation of cloud-cond ensation nuclei when the CR flux is high 
dSvensmark et al. 2007). Thus, an anti-correlation between tem- 
perature and CR flux i s to be expecte d and indeed repor ted by 
IShaviv & Ve"izer1d2003l) . Most recentlv. lSvensmarkl (1201 2l) found 
further evidence of a connection of nearby supernovae and their 
CR output and life on Earth . More detai l s on t he CR-climate 
connection can be found in IScherer et al.l (|2006). Although we 
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do not state here that we a dhere to this view in a ll aspects (see, 
e.g., the critical remarks in lOverholt et al.l l2009). we think that 
it gives a further interesting motivation to study the galactic CR 
distribution, and especially its longitudinal structure, in greater 
detail. 

The aim of this investigation is to calculate galactic CR spec- 
tra at different positions along the Sun's orbit around the galactic 
center and to analyse the influence of anisotropic diffusion on the 
longitudinal cosmic ray distribution. We first present the under- 
lying propagation model and its relevant input, like the diffusion 
tensor and its connection to the galactic magnetic field, the three- 
dimensional source distribution of CR and its connection to the 
spiral-arm structure and supernova (SN) occurrence, and loss- 
processes in the ISM. We also introduce our numerical solution 
method to the CR transport equation based on stochastic differ- 
ential equations. Finally, the calculated CR spectra and orbital 
flux variations are discussed and conclusion are drawn. Some 
earlier results on this topic can also be found in Effenberger et alj 
(1201 lh . 



2. The propagation model 

The basic transport theo ry of CRs is describ e d in many co ntem- 
porary monog raphs, e.g.|Schlickeiser ( 2002). IS tanevl (120041) and 
dShalchill20 09). Recentlv JStrong et al.1 d2007l) have surveyed the 
theory and experimental tests for the propagation of cosmic rays 
in the Galaxy. The considerations are com monly based on th e 
following parabolic transport equation (e.g., Ptus kin et al.l2 006): 
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where N{r, p, t) = p 2 f{r, p, t) is the differential intensity of CRs 
and / their phase space density, which is assumed to be isotropic 
in momentum space. As usual, r, p, and t denote the location 
in space, momentum and time and we use a galactic cylindri- 
cal coordinate system [r, ip, z]- The source term Q includes pri- 
mary particle injection, which, in this study, is considered to be 
only by supernovae and their remnant shock features. The spa- 
tial diffusion, in general, should be described by a tensor, but 
in most applications to galactic propagation so far, it is simpli- 
fied to a scalar coefficient k s , i.e. k = (a-, ; ) = (6ijK s ) (see the 
discussion above). An ordered motion of the IS M can be taken 
into account via the convection veloc ity u (e.g.. iFichtner et al] 
1199 It iPtuskin etal.lll997t IV61kll2007l) . but is neglected for this 
study due to its decreasing importance for higher CR energies. 
Continuous momentum losses are described by the momentum 
change rate p. Catastrophic loss processes like, e.g., spallation 
do not apply, since in this study only galactic protons are con- 
sidered. 



2.1. The anisotropic diffusion tensor 

As discussed in the introduction, generally, the diffusion of CR 
in magnetic fields with a prominent ordered field component is 
anisotropic with respect to this field orientation, i.e. stronger in 
field-parallel direction and weaker in the perpendicular direc- 
tions. This effect can be included in the propagation model by a 
diffusion tensor which is locally, that is, in a field-aligned coor- 
dinate system, diagonal: 

















Here, drift effects or aspects of non-axisymmetric turbulence 
(IWeinhorst et alJ |2008 ), which could lead to off-diagonal ele- 
ments in the diffusion tensor, are neglected. 

Since the CR transport is described in a global frame of refer- 
ence (i.e. the galactic frame with a cylindrical coordinate system 
in case of this study) the field-aligned tensor has to be trans- 
formed to this frame by the usual transformation 



Ak h A l 



(3) 



This transformation is analogous to the Euler angle transforma- 
tion known from classical mechanics. The matrix A = Rt,R2R\ 
describes three consecutive rotations Ri with A~ l = A T (since 
A € SO3). These rotations are defined by the relative orienta- 
tion of the local and the global coordinate system with respect to 
each other. 

If the two perpendicular diffusion coefficients are not equal, 
this transformation is of particular importance in establishing the 
appropriate orien tation in the calcu l ation o f the global diffusion 
tensor. Recently, Effenberger et al.l (I2012I) established a gener- 
alized scheme based on the local field geometry to account for 
this. In the present study, however, both perpendicular diffusion 
coefficients are set equal to reduce the set of unknown param- 
eters (connected to the unknown detailed turbulence properties 
in the ISM), i.e. k ± \ = k ± 2 = Furthermore, since the galac- 
tic magnetic field in consideration is to first order parallel to the 
galactic disc (see the discussion in the following subsection), the 
field tangential e, and the z-axis e z unit vectors provide, together 
with the completing third unit vector e„ = e z x e,, a well-defined 
coordinate system. These unit vectors represent the columns of 
the transformation matrix A. 

In order to complete the model of the diffusive part of CR 
propagation, the local tensor elements, i.e. k\\ and k ± have to 
be defined. For the parallel diffusion coefficient k\\, we assume 
the same broken power law depe ndence as has been taken fo r 
the scalar diffusion coefficient in Biisching & Potgieter (2008), 
namely: 



, P 



(4) 



(2) 



with a - 0.6 for/? > po, a = -0.48 for p < po, kq = 
0.027 kpc 2 /Myr and p = 4 GeV/c. 

Originally, this particular break energy of 3-4 GeV was 
motivated to fit the p lain diffusion model results like in 
Mos kalenko et al. (2002) to the obse rved Boron to Carb on ra- 
tios. A first refinement can be found in Ptuski n et al. (2006) and a 
more r igorous study has been performed by Shalchi & Biisching 
(l2010h . They confirm, by inclusion of turbulence dissipation ef- 
fects and the replacement of the quasilinear transport theory with 
a second-order diffusion theory, the possible existence of such a 
turnover in the parallel diffusion coefficient. 

The perpendicular diffusion k± is scaled to be a fraction of 
% i-e. 

K± = 6K|] (5) 

where the diffusion-anisotropy e is assumed to be i n the range of 
0. 1 to 0.01 for galactic protons with GeV energies (iShalchi et al.l 
2010). The actual variation of anisotropy with energy is an in- 
teresting aspect, but for the relatively small energy range up 
to 1 TeV, considered in this study, the anisotropy can be re- 
garded to first order as e nergy independent, thus following 
iGiacalone & Jokipiil d!999h . Since the latter author's finding of 
energy independence i s valid for both is otropic (see also the par- 
tially similar result bv lCasse et al.ll2002h and composite turbu- 
lence and given the to some extent different results found by 
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Shal chi et al.l (2010) assuming a Goldreich-Sridhar turbulence, 
conclusions on the detailed underlying turbulence properties 
should not be drawn on the basis of our study. 

2.2. Galactic magnetic field models 

As soon as anisotropic diffusion is considered, the knowledge of 
the large-scale magnetic field in the galaxy becomes importan t. 
Reviews on thi s su bject were written, e.g . , by Beck et alT dl996). 
iFerrierd (l200lh andHeiles & Haverkornl d2012l) . Pulsar rotation 



measure data (IHan e t al. 2006) give evidence for a counterclock- 
wise field orientation (viewed from the north Galactic pole) in 
the spiral arms interior to the Sun's orbit and weaker evidence 
for a counte rclockwise field in t he Perseus arm, see, however the 
criticism byfWielebinski (2005). In inter-arm regions, including 
the solar neighbourhood, the data suggests that the field is clock- 
wise. |HI3(l2QQ6j) proposed that the galactic magnetic field in the 
disk has a bi-symmetric structure with reversals on the bound- 
aries of the spiral arms. Magnetic fields in the general class of 
spira l galaxies were studied b y, e.g. JWielebinski & Beck! 12005) 
and iDettmar & Soidal (120061) and can be compared with that of 
our own galaxy. 

Fortunately, as long as drift effects are neglected, the actual 
sign-dependent orientation is of no relevance for the diffusion 
along and perpendicular to the magnetic field, which enables us 
to employ a simplified model of the galactic magnetic field that 
has no field reversals and is aligned to the spiral arm structure 
in the disc. Neglecting furthermore its weak halo-component, a 
simple model for the mean galactic magnetic field in cylindrical 
coordinates is given by 



B = Z?o(sin iff e r + cos iff e v ) - exp 



2cr 2 



(6) 



which is divergence free by construction. Here, if/ is the counter- 
clock logarithmic spiral arm pitch-angle, which has an a pprox- 
imate value of iff - 12°, according to the meta-study by Vallee 
J2001 . The same spiral arm parametrization is employed for the 
source distribution function discussed in the following section. 
Note, that the halo-scale parameter 0%,,, is not relevant in this 
context, since only the magnetic field direction is used for the 
construction of the diffusion tensor. 



2.3. Source distribution 

For the injection of CRs we assume a source distribution which 
follows the galactic spiral arm structure, where supposedly most 
of the supernovae (SN) occur. As a basis we take the same spi- 
ral arm model as m entioned above, i.e. the model established by 
Vallel<|2002l|2005), which consists of four logarithmic and sym- 
metrically positioned arms. Around these, we take a Gaussian 
shape (analogous to the approach in lShavivl d2003l) ) to yield an 
analytic expression for the source term Q, by summing up over 
all four arms (n € {1,2,3,4}): 




Fig. 1. Orientation of the galactic spiral arms in the present 
model. Norma, Scutum, Saggitarius, and Perseus arms are col- 
ored by green, blue, red, and purple, respectively. The black line 
shows the solar orbit and the galactic center region is marked in 
black as well. The four x-markings indicate the positions at 90°, 
108°, 126°, and 144°, where the CR spectra have been calculated 
(see Section[3]l. 




q n = Qop s exp 



(r-r„) 2 
2o-? 



2a-?/ 



(7) 



Fig. 2. Volume rendering visualisation of the input source distri- 
bution of CRs. The coloring indicates the source strength from 
red (low) to yellow (high) in arbitrary units. 



with r n = ro exp(k((p + </j„)). <p n introduces the symmetric rota- 
tion of each arm by 90°, i.e. ip„ = (n - \)n/2. k = cosiff with 
iff = 12° is the constant pitch-angle cosine of the spiral arms. 
Figure Q] illustrates the orientation of the spiral arms relative to 
the Sun's position and orbit. We take cr r - cr z - 0.2 kpc to 



have a reasonable inter-arm separation, while ro = 2.52 kpc ac- 
cording to Vallee's model. The model galaxy has the often as- 
sumed CYjmdric^lshape_with a radius of 15 kpc and a height of 
4 kpc (Buschin g et al.ll20 05) and the Sun's orbit is at a radius of 
r0 =7.9 kpc. The average spectral index s of the sources' power 
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Fig. 3. Sample paths of pseudo-particles in the galactic magnetic spiral field, coloured by three different colours for three particles, 
each starting at the same point in phase space (i.e. Earth's position at 1 GeV) and projected onto the galactic plane. The black lines 
show integrated magnetic field lines to illustrate the magnetic field orientation. The left panel shows sample paths for isotropic 
diffusion, with no visible effect of the magnetic field orientation. The right panel illustrates the preferential diffusion along the 
magnetic field for a simulation with anisotropy e = 0.1. Note that the exit point of the red particle in the right panel is actually the 
radial boundary, while the other particles all exit through the halo's z-boundary (not visible). 



law injection in momentum is set to s — 2.3 i n agreement with 
recent estimates on CR source spectra (see e.g. lPutze et alJ20l"il 
lAveetaU feOCm Figure [2] gives a visualisation of this source 
distribution. The overall source strength Qq is a free parameter 
which can be fitted to a given reference like a local interstellar 
spectrum (see also the discussion in Sectional. 

2.4. Loss processes 

The two most dominant loss processes for CR protons dur- 
ing their propagation through the ISM are energy losses due 
to pion-production for relativistic energies and ionization pro- 
cesses in the ISM plasma for lower energies (see Fig.l in 

iMannheim & Schlickeiserll 19941) 

According to Chapter 5 in [Schlickeiser (2002) the pion 
losses can be approximated for Lorentz factors y » 1 as 



- =1.4 



10- 16 («H/ 



2n H2 )A- 0A1 y L2 h- 1 



(8) 



where we assume a z-dependent ISM gas density with («/// + 



2n Hl ) = 



1.24 



cosh(30.'l 



(in units of particles per cm 3 and z in kpc, 



Busching & Potgieter 2008). The mass number A of a proton is 
just unity. A similar formula for the ionization losses is given by 



10 



1 Z 2 n„ 



P 



-P 3 



-eVc 



-i „-i 



(9) 



where the electron density is the same as the above gas density 
and the charge state of protons Z is equal to 1 . For the purpose of 
this study, the velocity factor/? = v/c of the particles is always 
much greater than the thermal electron B e which is related to x m 
by x m = (3/47r 1/2 ) 1/3 j8 c = 1.10/3 e . This means that the momen- 
tum loss rate scales approximately as p~ 2 . The total momentum 
loss rate entering Equation[T]is the sum of both loss processes. 



2.5. The numerical solution method based on stochastic 
differential equations 

To solve the transport equation for the problem setup introduced 
in this study, we employ a numerical solution scheme based on 
the It5 equivalence of a Fokker-Planck type equation and cor- 
responding stochastic differential equations (SDEs) involving a 
random Wiener process. This method has become increasingly 
popular in CR transport studies because of its numerical simplic- 
ity and conformance with modern computer architecture, i.e. its 
straightforward parallelization and scalability. Mentioning only 
a few examples, a starting point f or heliosphe ric studies of this 
kind can be found in the paper by Zhang! (1 1999b where he a pplied 
the method to CR mod ulation. More recently. |Pei et al.l (2010) 
and lStrauss et al.l d201 lb applied SDEs i n a more comprehensive 
heliospheric model. Farahat et aTl (2008) applied SDE s for a CR 



propag ation study in the Galaxy and, e . g.. iMarcowith & Kirkl 
(119991) as well as lAchterberg & Schurd ( 1201 lb calculated the 
shock acceleration of energetic particles. 

The basic idea in SDE schemes is to trace pseudo-particle 
trajectories from their origin forward in time or, alternatively, 
integrate backwards in time from the phase space point of inter- 
est. The particle trajectory is given by the integral of an SDE of 
the following form 



dxi — Ai(xi)ds + ^ Bij(xi)dWj 



(10) 



where the relation BB T = 0.5 k has to be fulfilled, that is, a root 
for the diffusion tensor k has to be determined. Here, dWj is a 
(multidimensional) Wiener process increment, which has a time- 
stationary normal-distributed probability density with expecta- 
tion value and variance 1. The deterministic part is directly 
related to the convection velocity in the transport equation, i.e. 
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Fig. 4. Calculated CR proton spectra (multiplied by E 2 ) at four 
different positions along the Sun's orbit (see Fig. Q] for the re- 
spective locations) in the galactic plane (z = 0), where 90° cor- 
responds approximately to the current solar system position and 
126° is inside of the Saggitarius arm. The upper panel shows the 
spectra for isotropic diffusion, while the middle and lower panel 
include anisotropic diffusion with e = 0.1 and e - 0.01, respec- 
tively. The latter spectra have been rescaled by factors of 0.15 
and 0.05 to account for the higher overall flux due to the stronger 
confi nement in the disk i n con trast to the isotropic case. The LIS 
from Web ber & Hig bie (2009) is plotted for comparison (black 
crosses). 



Aj = — m j. Numerically, this SDE is integrated via a simple Euler- 
forward scheme and the Wiener-process is simulated with the 



Box-Muller method (e.g.. [Box &Mullerlll958l) bv using 

dWi(s) = rj(s)^/ds (11) 

where rj(s) is equivalent to a Gaussian distribution N(Q, 1). 
The necessary random numbers are generated with the 
MTI 19937 version of the so called Mersenne Twister 
(iMatsumoto & Nishimural [l998l) . The integration parameter s is 
related to physical time by 



to 



(12) 



where fn is the final time for the backward method. The source 
contribution to the individual particle trajectory is added up by 
a path amplitude. Finally, in case of the backward method, all 
particle trajectories are weighted together to yield the resulting 
phase space density (i.e. the solution to the associated Fokker- 
Planck equation) at the starting phase space point. The boundary 
and initial conditions can be accounted for in the weighting, but 
for this study they are simply zero (corresponding to an escaping 
boundary condition for the CRs). We only apply the backward 
method in this study, since it is well-suited for the given problem. 
For more details on the numerical scheme and especially on the 
determin ation of the root o f th e diffusion tensor, w e refer the 
reader to lKopp et ail (1201 2|) and lStrauss et al.1 (1201 ll) where the 
basis of the code used in this study is discussed in greater detail. 

Exemplary pseudo-particle trajectories are shown in Fig. [3] 
There, the additional information contained in SDE calculations 
becomes obvious. The pseudo-particles' paths follow the field 
lines during their stochastic motion as soon as anisotropic dif- 
fusion becomes relevant. Consequently, the modification to the 
diffusion process becomes directly visible in such trajectories. 
However, one has to keep in mind that these are not real particle 
trajectories or gyro-center motions, but only tracers of the phase 
space of the diffusion-convection problem. 

3. The resulting spatial and spectral CR distribution 

We calculated CR proton spectra within the introduced model at 
four different positions along the Sun's galactic orbit. The posi- 
tions are indicated in Figure Q] We took a very long integration 
time (fo ~ 10000 Myrs) to assure that we approached a steady 
state situation, which is confirmed by checking that all particles 
have exited the computational domain. For each phase-space 
point, 10 4 pseudo particle trajectories have been computed. A 
comparison between the calculated spectra in the case of pure 
isotropic diffusion (upper panel) and two anisotropic case, with 
weak (e = 0.1, middle panel) and strong (e = 0.01, lower panel) 
diffusion anisotropy, is shown in Figure [4] For comparison, 
the lo cal interstellar spectrum (LIS) given by Webber & Higbie 
(2009) is included in the plots (WH0 9 he r eafter) . We have used 
the parametrization given in iHerbst et al. l (f20T0h . where a com- 
parison between a few proposed LIS can be found as well. In 
face of the still imprecisely known modulation effects on mea- 
sured spectra inside the heliosphere (see, e.g., iFlorinski et alj 
l201lHScherer et alJ201ll) . such an LIS parametrization can give 
only a rough orientation on what to expect for galactic CR prop- 
agation studies. Our results have been rescaled to fit approxi- 
mately to the WH09 LIS in the isotropic diffusion case, by ac- 
counting for the free parameter Qq in the source strength. The 
anisotropic spectra have been rescaled again, respectively. To 
yield the good agreement shown in both upper panels of Figure|4] 
between the calculated spectra and the WH09 LIS, the break in 
the diffusion coefficient introduced above as well as both con- 
tinuous loss processes are required. The inclusion of the latter is 
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an im provement over earlier studies, like [Biisching & Potgieter 
(2008), where only a parametrized catastrophic loss term was 
considered. 

The spectra for different positions along the Sun's galactic 
orbit show only very little variation in the isotropic diffusion 
case. Particularly, the variation is largely independent of energy 
over the entire energy range considered. In contrast to this, the 
variation is much stronger for the anisotropic case, depending on 
the imposed diffusion anisotropy. The differences are, in these 
cases, dependent on energy as well. For high energies, the spec- 
tra start to converge again towards the isotropic differences. This 
is due to the increasing dominance of escape losses for these 
high energies. For lower energies, the pion and ionization losses 
are much more important than in the isotropic case, because the 
confinement time of CRs is longer as a result of the reduced dif- 
fusion perpendicular to the disk. 

Notably, the spectrum at Earth for the weak anisotropic case 
fits even better to the reference LIS than the pure isotropic re- 
sult, which shows that, depending on the overall parameter set, 
anisotropic diffusion scenarios can improve on the model results 
of conventional studies. In this context, one has to keep in mind 
however, that the precise form of the low energy LIS and its 
connection to the galactic spectrum on a kpc scale is yet unclear 
and depends on modulation effects in the heliosphere as well as 
similar effects i n the local solar sy stem environment (see, e.g., 
the discussion in lScherer et al.ll20lTl) . Furthermore, the assumed 
break in the diffusion coefficient may be different or even ab- 
sent in a more complete propagation scenario, since up till now 
it is mainly phenomenologically motivated, to yield the expected 
local spectra. 

The spectra for the strong anisotropic case (e = 0.01) show 
some significant deviations from the expected spectral shape due 
to the largely increased relative importance of the loss processes, 
resulting e.g. in a flatter high energy spectrum. In the context of 
the model setup of this study, this means that such a high dif- 
fusion anisotropy is probably unrealistic. Nevertheless, we in- 
cluded this case since it shows the resulting large orbital varia- 
tion at lower energies (see also Fig. [5) where the spectral shape 
is still unclear. In addition, models with different structures in 
the galactic halo, namely with different gas densities and halo 
heights, as well as a possible magnetic field component perpen- 
dicular to the disc, may alter the resultant spectra further, due to 
a changed influence of the loss processes. These aspects could 
be further clarified in a subsequent study which takes different 
CR species and more sophisticated magnetic field models into 
account. 

In Figure [5] the orbital flux variation along the Sun's or- 
bit is plotted against longitude, to illustrate further the varia- 
tion with longitudinal position. We consider two different ener- 
gies, namely 1 GeV and 100 GeV and two galactic distances of 
r — 5 kpc and tq =7.9 kpc, again for isotropic (left panel) and 
strong anisotropic (e = 0.01, right panel) diffusion. The inclu- 
sion of a second radius at only 5 kpc is motivated by the recent 
claim that the Sun may have migrated outw ards during its life- 
time in the galaxy (Nieva & Przvbilla 2012). It can be seen that 
the variation is much weaker in a closer galactic orbit and has 
a different phase, as a result of the smaller inter-arm separation. 
For all cases, the overall shape of the variation is not a simple 
sinusoidal profile due to the non-perpendicular transit of the Sun 
through the arms (see again Figure [TJ. The amplitude of varia- 
tion is much more pronounced for the anisotropic case, that is, it 
can be as large as a factor of 6, while in the isotropic case it is 
only a factor of about 2. An increased diffusion anisotropy e will 
increase this difference even further. 



4. Conclusions 

In this study the effects of anisotropic diffusion of galactic CR 
protons have been analyzed. For the solution of the steady-state 
diffusion equation a numerical method based on stochastic dif- 
ferential equations has been used, which also accounts for en- 
ergy loss processes. The computed spectra along the Sun's galac- 
tic orbit show larger variations for the anisotropic cases when 
compared to the scalar diffusion model. Furthermore, for the 
chosen parameters, a moderate diffusion anisotropy (e = k x /k\\ = 
0.1) leads to a result which, in our setup, is more consistent with 
recent estimates of the local interstellar proton spectrum than the 
results for purely isotropic diffusion. 

We therefore conclude that the diffusion tensor as well as the 
CR source distribution is an important feature in determining the 
solution of the transport equation in a three-dimensional model 
of galactic CR p r opa gation. This result fits well into the findings 
bv lHanasz et al.l (120091) claiming that anisotropic diffusion is an 
essential requirement for the CR driven galactic dynamo effect. 
Additionally, these results imply that i n the context of the re- 
cently proposed CR-climate connection (IShaviv & Veizeif 2003) 
the expected CR flux variation may be even larger than previ- 
ously estimated, although at present, this topic is still highly 
speculative. 

We emphasize that the introduced model indicates further 
research opportunities by including additional effects like a vari- 
able spectral source index and time-variable CR sources, de- 
pending on supernova type. The necessary time-dependent cal- 
culations are in principle possible with the current mode l setup, 
and the earlier work by Biisching & Pot gieterl (120081) shows 
some results on this already. Furthermore, the consideration of 
additional CR species and their relevant loss processes may yield 
further insight to constrain the transport parameters. Finally, this 
study points at the necessity to acquire more detailed models of 
the galactic magnetic field, to asses its impact on the transport 
processes of CRs. 
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